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Abstract 

The transformation stretch tensor plays an essential role in the evaluation of conditions of compatibility 
between phases and the use of the Cauchy-Born rule. This tensor is difficult to measure directly from 
experiment. We give an algorithm for the determination of the transformation stretch tensor from x-ray 
measurements of structure and lattice parameters. When evaluated on some traditional and emerging phase 
transformations the algorithm gives unexpected results. 
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The structural transformations commonly occur in application of functional materials. Typical 
examples of phase transformation driven phenomena include shape memory alloys, ferroelectric- 
ity, piezoelectricity, colossal magnetoresistance and superconductivity. It has been demonstrated 
that material reliability depends, essentially, on the reversibility of the transformation. It is there¬ 
fore important to understand how reversibility can be achieved and how the transformation occurs 
at the lattice and atomic level. The transformation stretch tensor , U, is the stretch part of the 
linear transformation that maps the crystal structure from its initial phase to the final phase [1-4]. 
Recently, the reversibility, the thermal hysteresis, and the resistance to cyclic degradation of func¬ 
tional materials have been linked to properties of the transformation stretch tensor. For example, 
when the middle eigenvalue Xj of U is tuned to the value 1 by compositional changes, the mea¬ 
sured width of the thermal hysteresis loop drops precipitously to near 0 in diverse alloy systems 
[5-7]. Assuming the Cauchy-Bom rule for martensitic materials [3, 8, 9], the condition X 2 = 1 
implies a special condition of compatibility between phases by which the undistorted austenite 
phase and a single undistorted variant of the martensite phase fit perfectly together at an interface. 
Even stronger conditions of compatibility known as the cofactor conditions (Xi — 1 together with 
either |U _ 1 e| = 1 or |Ue| = 1, where e is unit vector on a 2-fold symmetry axis of austenite), lead 
to even lower hysteresis and significantly enhanced reversibility during cyclic transformation [10]. 
U also plays an important role in determining the elastically favored orientations of precipitates 
for diffusional transformations [4, 11]. 



FIG. 1. Non-uniqueness of Cauchy-Born deformation gradient from (a) square lattice to (b) oblique lattice. 
Red, blue and green balls represent different atomic species. Gray dots are lattice points 
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In principle, the determination of the stretch tensor U for a structural transformation is straight¬ 
forward. Suppose the primitive lattice vectors of initial and final phases are, respectively, linearly 
independent vectors a, and b, for i = 1,2, ...d where d is the dimension of the lattice. A nonsingular 
linear transformation F can be defined uniquely by 

Fa, = b,, i=l,2,...d, (1) 

and the polar decomposition of F is written F = QU, where Q is orthogonal and U is positive- 
definite and symmetric, called the transformation stretch tensor. The notation a, —» b, denotes 
the lattice correspondence. In the case of transformation in Fig. 1, one choice of the lattice 
correspondence can be ai —> bi, a 2 —> b 2 where ai = [1,0], a 2 = [0,1] and bj = [a,0] and b 2 = 
[Z? cos /3 , sin /3 ]. 

As is well-known [12, 13], F and U are not uniquely determined by the two lattices. This fol¬ 
lows from the fact that there are infinitely many choices of lattice correspondence. From Fig. 1, the 
alternative set of vectors ai and ai +a 2 describes the same lattice (a), which results in a different 
correspondence from (a) to (b). This obviously changes the F and thus the transformation stretch 
tensor U. More generally, any two sets of primitive lattice vectors for a given lattice are related by 
a lattice invariant transformation [2] i.e., a unimodular matrix of integers. If we allow an invariant 
transformation for both initial and final phases, the ambiguity of F is F —» A^FA^ 1 where A^ 
and A(f) denote the lattice invariant transformation for initial and final lattices, respectively. 

The linear transformation F represents the change of periodicity of the two phases. The indi¬ 
vidual atoms denoted by the red, blue and green balls in Fig. 1 may shuffle in various ways, giving 
rise to different space group symmetries, but it is the linear transformation F that relates to macro¬ 
scopic deformation and therefore to conditions of compatibility [1, 7-9, 12, 14-18]. This idea is 
formalized by the weak Cauchy-Born rule [8, 19]. This rule is used to define the dependence on 
deformation of the free energy at continuum scale from the free energy density at atomistic scale 
for complex lattices with multiple atoms per unit cell and inhomogeneous deformations. Inho¬ 
mogeneous deformations y(x) locally satisfy the same rule as above: b, = Vya,, where a,- and 
b, represent the local periodicity. Note that we use a geometrically exact description here. A 
geometrically linear description (i.e., as in linear elasticity) would not be sufficiently accurate to 
describe transformations here for the purposes of imposing the conditions of compatibility (see [7] 
for calculations of the error in various cases). 

Based on a natural intuition that “a mode of atomic shift requires minimum motion” [20], Bain 
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proposed a famous lattice correspondence in 1924 for the formation of bee aFe from fee yFe . 
The correspondence has been well-accepted and applied to study numerous phase transformations 
[2, 3, 21-25]. To illustrate how easy the Bain correspondence misses the smallest strain, we 
construct an example of transformation from a bee lattice with ciq — 1 to a monoclinic lattice with 
a = 0.96l,b— 1.363 ,c— 1.541,and/3 = 97.78°. Fig. 2(a) shows the bee lattice with two sublattice 
unit cells (red and blue). Conventional wisdom would say that the Bain correspondence (red —>■ 
gray in Fig. 2(b), bottom) is appropriate for this transformation. However, our algorithm proposed 
later in this letter reveals an unexpected alternative correspondence (blue —>• gray, Fig. 2 (b), top). 
Both contain 4 lattice points (n — 4) in the unit cells, and the shape and size of them are similar 
to the primitive cell of monoclinic lattice. Fig. 2(b) shows the comparison of distortions for both 
transformation mechanisms. Notice that both mechanisms give exactly the same final monoclinic 
lattice. However, by quantitative calculation, the principle strains for the new correspondence are 
in fact smaller than those for the Bain correspondence. 



FIG. 2. The least atomic movements during the structural transformation, (a) The bee lattice and two of its 
sublattices (red and blue) of size 4. (b) Comparison between these bee sublattice unit cells and the primitive 
cell of the final phase (gray; for clarity atoms in the unit cell are not shown). 

The significance of finding the correct lattice correspondence for structural phase transforma¬ 
tions is emphasized in the literature [12, 13]. The problem was well-appreciated by Lomer [26] 
as early as the mid-1950s. In his study of the mechanism of the j8 —> a phase transformation of 


4 



















































U 9 8.6Cr L4 , he examined theoretically (by hand) 1,600 possible transformation mechanisms, and 
reduced this to three correspondences having the smallest principle strains, which he considered 
the likely candidates. 

Direct experimental measurement of the macroscopic finite strain of transformation, together 
with accurate structural characterization by X-ray diffraction provides a possible way to determine 
the lattice correspondence and thus the transformation stretch tensor. But this is technically dif¬ 
ficult due to (i) the need for an oriented single crystal, (ii) the need to remove the inevitable fine 
microstructures that form during transformation due to constraints of compatibility, and (iii) the 
need for an accurate measure of full finite strain tensor along known crystallographic directions. 
We also noticed that using a state-of-art high resolution TEM on a pre-oriented single crystal sam¬ 
ple can not definitively remove the ambiguities among many lattice correspondences due to some 
inevitable obstacles: tracking the evolution of diffraction spots in a fast structural transformation 
process, simultaneously indexing both phases, and most significantly, finding a special zone that 
can unambiguously reveal the differences among various lattice correspondences. 

In this letter we propose an algorithmic approach to search the N best choices of lattice cor¬ 
respondence for a structural transformation, by minimizing a particular strain measure between 
initial and final lattices. The input to the algorithm is the underlying periodicities (the remaining 
space group information is not needed) and the lattice constants of the two phases. The output 
from the algorithm is the N best choices of lattice correspondence and the associated transforma¬ 
tion stretch tensors. Users can customize how many solutions they like by manipulating N. The 
results can be used as a reference by the advanced structural characterization facilities for the de¬ 
termination of orientation relationships, and it can be integrated with first principles calculations 
to give starting points for the determination of energy barriers or interfacial distortion profiles. 

Consider a Bravais lattice ££ = {Xu'e, : n 1 ,... n d G Z d } determined by linearly independent 
lattice vectors ei,... ,e^ G W 1 , i — 1,... ,d, and assemble the lattice vectors as the columns of a 
dxd matrix E = (ei,... ,e</). ££ can equivalently be denoted 

& = Jf(E) = {r G R d : r = E£,£ G Z d ). 

Without loss of generality, by switching the sign of ei if necessary, we assume that detE > 0. This 
determinant is the (J-dimensional) volume of a unit cell of J§?(E). 

Given two lattices Jf(E) and Jf(E'), the dxd nonsingular matrix L satisfying E' = EL is 
called the correspondence matrix from jSf(E) to jSf(E'). As noted above, the two lattices Jf(E) 
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and Jf(E') are the same if and only if the correspondence matrix L is a unimodular matrix of 
integers, or, briefly, L 6 GL(d. Z). If a correspondence matrix L is a matrix of integers with 
| detL| > 1, then J?f (E r ) is a sublattice of Jf(E). The quantity | detL| is the volume ratio of the 
unit cell of Jf (E') to that of Jf (E). 

Correspondence matrices are often reported for conventional rather than primitive descriptions, 
particularly for 7 of the 14 types of Bravais lattices in 3D. For example, the conventional descrip¬ 
tion for an fee lattice with lattice parameter ao is an orthogonal basis, so E con v = «()I = Ex, where, 
for example, 


1 0 1 


1 1 -1 

1 1 0 

, x = 

-1 1 1 

0 1 1 


1 -1 1 


Here, det% = 4 so the volume of the conventional unit cell is 4 times that of the primitive cell. 
From now on, the symbol x is reserved for a correspondence matrix from the primitive to conven¬ 
tional unit cell of a Bravais lattice: E conv = E%. 

We seek a sublattice of (E A ) that is mapped to the primitive lattice of Eg). (The algorithm 
can easily handle the case in which we take sublattices of both lattices.) As above, let E A = 
(ai, ...,a c i) and Eg = (bi, ...,b^). Let £ e Z dxd , det£ > 0, be the correspondence matrix giving 
the sublattice Jf(E^£) that is mapped to the final lattice Jf( Eg) during the transformation. The 
basic equation (1) in this case becomes FE A £ = Eg, and the transformation stretch tensor U is the 
unique positive-definite square root of F r F. 

We introduce the following function as a measure of the distance from U to I: 

dist(^E A ,Eg) = ||(F r F)- 1 -l|| 2 

= \\e a £e~ b 1 e- b t £ T E T A -i\\f. ^ 

|| ■ || denotes the Euclidean norm, ||A|| = Vtr A 7 A. The distance (2) is independent of rigid rota¬ 
tions of both lattices, and is particularly attractive from the point of view of symmetry. Physically, 
it represents the Lagrangian strain of the structural transformation. The use of inverse of F y F 
avoids possible noninvertibility of £ that may arise during the minimization process. In addition, 
this norm is exactly preserved by point group transformations of both Bravais lattices. That is, 
if orthogonal tensors R A and Rg are, respectively, in the point groups of 2§f(E A ) and jSf(Eg), 
i.e., 2£?(E a ) = 22 ?(R a E a ) and J§?(Eg) = =£?(RgEg), which, by the above implies that there ex¬ 
ist associated matrices /r A an d /rg such that R A E A = E A /r A and RgEg = Eg/rg then the distance 
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transforms as 


dist (haZhb, Ea , Eg) = dist(4 E A , Eg). (3) 

Note that /i A ,g are integral matrices of determinant ±1, so detf! = dct /r A ^/rg. Thus, immediately 
one minimizer of the distance with assigned determinant gives the expected symmetry-related min¬ 
imizes. Physically, in the typical case of a symmetry-lowering transformation, e.g. the marten¬ 
sitic transformation, the distance function (2) automatically gives the equi-minimizing variants of 
martensite. 

As noted above it is typical to report the correspondence matrix in terms of the conventional 
basis instead of the primitive one. If £.* is a minimizer of dist(CE A .Eg) the conversion is done by 
L* = x2 l FXB. N° te that L* is not necessarily a matrix of integers. 

A significant property of the distance function (2) will be used to justify our algorithm below. 
Fixing E a and Eg, the distance function can be trivially extended to a function over real matri¬ 
ces, /(L) = dist(L,E A ,Eg). Denoting Xg = E A LE^ 1 E B r L r Ej and using Xg ■ I < ||Xg|| ||I|| = 
y/3 ||Xg||, we have 

/( L) = ||X l || 2 —2X L -I + 3 

(4) 

^||X i || 2 -2v / 3||X t ||+3 = (||X I .||-v / 3) 2 , 

Choose any integral matrix i\ and define Ci = f{t\ ). By (4) the minimizer(s) of /( L) necessarily 
lie in the bounded set ||Xg|| < a/3 + \/Ci, that is, ||Xg|| 2 < 3 + Ci + 2\/3Ci. Let a be the minimum 
of ||Xg|| 2 under the constraint ||L|| = 1, then we have 

a||L|| 4 < ||X L || 2 < 3 + Ci +2 a/3C)". (5) 

That is, all the L’s such that /(L) < C live in the sphere with the radius of ((3 + Ci -t-2 v / 3C7)/a) ! / 4 
in M 9 . 

Here is a brief outline of the algorithm for the determination of the N best transformation stretch 
tensors and their associated lattice correspondences: 

1. Calculate the primitive bases and the transformation matrices for the conventional cells from 
the input lattice parameters: Ea, b and xa, b- Calculate a by minimizing the term Xg with 
respect to L for all ||L|| = 1. 

2. Choose N integral matrices iu i = 1,..., N as the initial guess of the solution list such that 
det £j is close to detEg/dctE A and dist(T,-,E A ,Eg) is small. 
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(a) 




FIG. 3. Two possible lattice correspondences in an FCC to monoclinic transformation, (a) (010) projection 
of the FCC lattice: the dark (resp. light) atoms are in the y = 0 (resp. y = 1/2) planes. The solid blue 
and red lines represent the two lattice correspondences respectively for m = 4, where the the Bain corre¬ 
spondence is in blue. The dashed blue lines indicate the modulation numbers m = 1,2,3,4. (b) shows the 
dependence of the values of the distance function on the modulation of the monoclinic c-axis for the two 
lattice correspondences. 

3. Let Ci be the maximum /(•£,) for s in the solution list. 

4. Calculate the distance for all integral matrices in the sphere of radius of ((3 + Ci + 
2-v / 3Ci)/a) 1 / 4 . Update the solution list as necessary. If the solution list is changed, re¬ 
peat from step 3. 

5. For each solution £j, calculate the Cauchy-Born deformation gradient F, = E^E^,-) -1 and 
the transformation stretch tensor U; = (F^F,-) 1 / 2 . Finally, rewrite all the solutions in the 
conventional bases: L* — Xa^X B- 

Note that the algorithm converges in a finite number of steps and gets all matrices with the N 
lowest distances (up to the degeneracy in (3)) because it searches through all matrices of integers 
satisfying the rigorous bounds (5). 

In Fig. 3 we give an example computed by the algorithm that reveals a switch from Bain corre¬ 
spondence to a new correspondence with increasing lattice complexity. Consider a transformation 
from an fee lattice with lattice parameter ao = 2 to a monoclinic lattice with lattice parameters 
a — 1.41, b — 1.99, c = 1.42 m, /3 = 86°, where the integer m > 0 denotes the modulation along 
monoclinic c-axis. Fig. 3(a) shows the undeformed fee lattice projected onto (010) plane. The 
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two correspondences given by the algorithm are depicted for the m = 4 case. Fig. 3(b) shows the 
change in distance function for the two correspondences with m varying from 1 to 16. Initially 
Bain correspondence is much smaller than the new one, however it loses its privilege after the 
7th modulation. The results suggest that both kinds of lattice correspondence can be feasible in a 
structural transformation for some special lattice parameters, and in this case m = 7 has this spe¬ 
cial status. As mentioned above, these long stacking period structures are common in martensitic 
phase transformations. 

Table I shows the results calculated by the algorithm for six materials. The types of transfor¬ 
mation are diverse and the principle stretches are consistent with the references. Among these 
examples, we list two solutions for Zn 45 AU 30 CU 25 . The material has been recently found to sat¬ 
isfy the cofactor conditions (the 2 constraints on U explained in paragraph 1) [7], which have 
been shown [10] to promote unusually low thermal hysteresis (« 2°C) and enhanced reversibil¬ 
ity, owing to a fluid-like flexible martensite microstructure. It was believed [10] to transform 
by the second solution, Table I. However, the first solution is the one having the smallest trans¬ 
formation strain. Coincidentally, the new transformation stretch tensor also satisfies closely the 
cofactor conditions. To investigate this further, the same sample of Zri/^ AU 30 CU 25 used in [10] 
was characterized by synchrotron X-ray Laue microdiffraction. The experiment has been con¬ 
ducted on beamline 12.3.2 of the Advanced Light Source, Lawrence Berkeley National Labora¬ 
tory. Details on the experimental setup can be found in [27]. The Laue patterns were collected 
continuously as heating/cooling through the transformation temperature. These patterns were an¬ 
alyzed and indexed using the XMAS software [28]. The orientation relationships are determined 
as the closest parallelisms of the crystallographic planes and zone axes between the indexed Laue 
patterns of austenite and martensite respectively. They are (206) a ||(2034) m , (204) a ||(1026) m , 
[211] a | | [269 l] m , [010] a 11 [010] m and [ 110] a 11 [891] m (see supplementary for indexed diffraction pat¬ 
terns). However, this determination with accepted error bars does not definitively distinguish these 
two mechanisms, since these relationships are so close that one could imagine that both mecha¬ 
nisms occur simultaneously in the material. 

In addition to the reversible martensitic transformation, the algorithm is applicable to a wide 
range of phase transformations even if the initial and final crystal structures do not have a 
group/sub-group relation. Examples are Tic^Mn^ and Sb 2 Te 3 /PbTe (Table I). The algorithm can 
be also applied to organic materials when the molecular chains have sufficient periodicity. One 
extreme example is the polymorphic transformation between two triclinic lattices in terephthalic 
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TABLE I. Transformation principle stretches (p. s.), the associated lattice correspondences (lat. cor.) and 


derived orientation relationships (o. r.) for various phase-transforming materials 


materials 

p. s. 

lat. cor. 

derived 0 . r. 


0.9363 

[j0j] L2l ->• [100] M 

(204) L2l ||(1 026 ) m 


1.0017 

s 

0 

S 

t 

rT 

O 

s 

[1 10] L2i 11 [891] m 

Zn45 AU 30 C 1125 [ 10 ] 

1.0589 

[405] L2l ->• [001] M 

[21 l] L2l 11 [2691] M 

L2j -5>M18R 

0.9363 

[j 0 j] L2i -»■ [100]m 

(204) L2i 11(1 027) m 


1.0006 

[010] L2l -> [010]m 

[110]l2j [991] m 


1.0600 

[|o|] L2l -> [001]m 

[21 l] L2l 11 [27 91] M 


0.9178 

[|0|]a^[100]b 

(H0)/3j 11(121)/ 

CuAboNL [29] 





1.0231 

CO 

0 

0 

t 

c 

0 

0 

[Ullft II[210]/ 

Pi ->/ 





1.0619 

[jOjJa — >■ [001] B 



0.9052 

[010] c -»■ [100]h 

(111 )c 11 (214) h 

Ti 95 Mn 5 [16] 





1.0164 

[iii]c^[010] h 

[121] c ||[201] h 

bcc —> hexagonal 





1.1086 

[101] c [001] h 



0.9791 

[11% ->[100]/»* 

(100)^11(111)^ 

RusoNbjo [30] 





1.0024 

[110]p/->[010] r 

[01 lift II [110]^ 

P'^P' 





1.0169 

[11%->[00%, 



0.9384 

[jjOJc —t [100] h 

( 110 )c (oio) h 

Sb 2 Te 3 / PbTe [4] 





0.9384 

[oH] c -> [oio] h 

[001]e| [48l] h 

fee —> hexagonal 





1.0779 

[222] c — > [001]h 



0.8244 

[012]! [ 100]n 

[100]i||[112]n 

Terephthalic acid [31] 





0.9373 

[110]i [010]n 

[010]i| | [102]n 

triclinic I —> triclinic II 





1.3424 

[001 ] i [001] 11 

[110]i||[010]n 


acid (see Table I). In this case the calculated principle stretches agree well with the measured 
macroscopic deformation of the polymorphic transformation of this material. 
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